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Abstract. We investigate nonlocal interaction equations with repulsive-attractive radial 
potentials. Such equations describe the evolution of a continuum density of particles in 
^ ^ which they repulse each other in the short range and attract each other in the long range. 

We prove that under some conditions on the potential, radially symmetric solutions converge 
exponentially fast in some transport distance toward a spherical shell stationary state. Oth- 
erwise we prove that it is not possible for a radially symmetric solution to converge weakly 
toward the spherical shell stationary state. We also investigate under which condition it is 
possible for a non-radially symmetric solution to converge toward a singular stationary state 
I I supported on a general hypersurface. Finally we provide a detailed analysis of the specific 

case of the repulsive-attractive power law potential as well as numerical results. 

< 

^ 1. Introduction 

> 

a Nonlocal interaction equations are continuum models for large systems of particles where 
, every single particle can interact not only with its immediate neighbors but also with particles 

far away. These equations have a wide range of applications. In biology they are used to model 
the collective behavior of a large number of individuals, such as a swarm of insects, a flock of 
^ birds, a school of fish or a colony of bacteria [MllMlEZlESlETlEZllMlEailllESlE 

lO [151 [m [13]. Ill these models individuals sense each other at a distance, either directly by sound, 

CN sight or smell, or indirectly via chemicals, vibrations, or other signals. Nonlocal interaction 

equations also arise in various contexts in physics. They are used in models describing the 
0^ evolution of vortex densities in superconductors [6ll[Hl[5lll2l[2l[Illl3l[3ll[ll]. They 

also appear in the modeling of dynamics of agglomerating particles in two dimensions (with 
I loose links to the one-dimensional sticky particles system) [l8]. They also appear in simplified 

inelastic interaction models for granular media [28[ \20\ \58[ HT] . Going back to biology, nonlocal 
• '-j interaction equations arise also in the modeling of the orientational distribution of F-actin 

/\ filaments in cells [351 [381 [5l| . 

In their simplest form, nonlocal interaction equations can be written as 

^ + div{fiv) = , v = -VW * n (1) 

where x) = ^t{x) is the probability or mass density of particles at time t and at location 
X E M^, W : — )• M is the interaction potential and f (t, x) is the velocity of the particles. 
We will always assume that the interaction potential W{x) = k{\x\) is radial and C^- or 
C^-smooth away from the origin, depending on the results. Typically the potentials we will 
consider have a singularity at the origin. 

When the potential W is purely attractive, i.e. is a radially symmetric increasing 
function, then the density of particles collapse on itself and converge to a Dirac Delta function 
located at the center of mass of the density. This Dirac Delta function is the unique stable 
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steady state and it is a global attractor [18j. The collapse toward the Dirac Delta function can 
take place in finite time if the interaction potential is singular enough at the origin and several 
works have been recently devoted to the understanding of these singular measure solutions 



In biological applications however, it is often the case that individuals attract each other in 
the long range in order to remain in a cohesive group, but repulse each other in the short range 
in order to avoid collision [47\ H9] . This lead to the choice of a radially symmetric potential 
W which is first decreasing then increasing as a function of the radius. We refer to these type 
of potentials as repulsive-attractive potentials. Compared with the purely attracting case 
where solutions always converge to a single Delta function, nonlocal interaction equations 
with repulsive-attractive potentials lead to solutions converging to possibly complex steady 
states. As such, nonlocal interaction equations with repulsive-attractive potentials can be 
considered as a minimal model for pattern formation in large groups of individuals. 

Whereas nonlocal interaction equations with purely attractive potential have been inten- 
sively studied there are still relatively few rigorous results about nonlocal interaction equa- 
tions with repulsive-attractive potential. The ID case has been studied in a series of works 
|33l [32l [52] . The authors have shown that the behavior of the solution depends highly on 
the regularity of the interaction potential: for regular interaction, the solution converges to a 
sum of Dirac masses, whereas for singular repulsive potential, the solution remains uniformly 
bounded. They also showed that combining a singular repulsive with a smooth attractive 
potential leads to integrable stationary states. Pattern formation in multi-dimensions have 
recently been studied in |39| I60j . In these two works, the authors perform a numerical study 
of the finite particle version of ([T]) and show that a repulsive-attractive potential can lead 
to the emergence of surprisingly complex patterns. To study these patterns they plug in 
([T]) an ansatz which is a distribution supported on a surface. This give rise to an evolution 
equation for the surface. They then perform a linear stability analysis around the uniform 
distribution on the sphere and derive simple conditions on the potential which classify the 
different instabilities. The various instability modes dictate toward which pattern the solu- 
tion will converge. They also check numerically that what is true for the surface evolution 
equation also holds for the continuum model ([T]). In another recent work jM] the specific case 
where the repulsive part of the potential is the Newtonian potential and the attractive part 
is polynomial is analyzed showing the existence of radially compactly supported integrable 
stationary states. They also study their nonlinear stability for particular cases. 

In this paper we focus primarily on proving rigorous results about the convergence of 
radially symmetric solutions toward spherical shell stationary states in multi-dimensions. 

Definition 1 (Spherical Shell). The spherical shell of radius R, denoted 6r, is the probability 
measure which is uniformly distributed on the sphere dB{0,R) = {x G M''^ : \x\ = R}. 

Given a repulsive-attractive radial potential whose attractive force does not decay too fast 
at infinity, there always exists an i? > so that the spherical shell of radius i? is a stationary 
state as it will be remarked below. One need then to address the question of wether or not 
this spherical shell is stable. It is classical, see [31 [201 EHl [21]) tliat the equation ([T| is a 
gradient flow of the interaction energy 
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with respect to the euclidean Wasserstein distance. Thus, stable steady states of ([T]) are 
expected to be local minimizers of the interaction energy. Simple energetic arguments will 
show that in order for the spherical shell of radius i? to be a local minimum of the interaction 
energy, it is necessary that the potential W satisfies: 

(CO) Repulsive-Attractive Balance: uj{R,R) = 0, 

(CI) Fattening Stability: diuj{R,R) < 0, 

(C2) Shifting Stability: diu{R,R) + d2Uj{R,R) < 0, 

where the function w : — ^ M is defined by 

uj{r,rj) = — —[ VW{rei - rjy) ■ eida{y), (2) 

(^N JdB{0,l) 

(Tat is the area of the unit ball in M^, ei is the first vector of the canonical basis of M^, 
da denotes the volume element of the manifold where the integral is performed and = 
(0, +00) X (0, +00). Condition (CO) simply guarantees that the spherical shell 6r is a critical 
point of the interaction energy. We will see that if condition (CI) is not satisfied then it is 
energetically favorable to split the spherical shell into two spherical shells. Heuristically this 
indicate that the density of particles, rather than remaining on the sphere, is going to expand 
and occupy a domain in of positive Lebesgue measure. If condition (CI) is not satisfied 
we will therefore say that the "fattening instability" holds. It can be easily checked that if 
uj{R, R) = 0, then diU}{R, R) is simply the value of the divergence of the velocity field on the 
sphere of radius R. So the fattening instability corresponds to an expanding velocity field 
on the support of the steady state. We will also see that if condition (C2) is not satisfied 
it is energetically favorable to increase or decrease the radius of the spherical shell. This 
instability will be referred as the "shift instability" . 

We now outline the structure of the paper and describe the main results. In the preliminary 
section, section [2| we derive (C0)-(C2) from an energetic point of view and we show that 
they correspond to avoiding the fattening and shift instability. We also study the regularity of 
the kernel uj defined by ([2]) . A good understanding of the regularity of a; will be necessary for 
later sections. We also remind the reader of previous results from [HlHj about well posedness 
of ([1]) in LP(M^). Section [3] is devoted to a detailed study of the fattening instability, both 
in the radially symmetric case and in the non-radially symmetric case. We first show that if 
condition (CI) is not satisfied then it is not possible for a radially symmetric L^-solution to 
converge weakly-* as measures toward a spherical shell stationary state. We then investigate 
singular stationary states supported on hypersurfaces which are not necessarily spheres. Such 
steady states have been observed in numerical simulations [391 I60j . We show that if the 
divergence of the velocity field generated by such stationary state is positive everywhere on 
their support, then it is not possible for an L^-solution to converge toward the stationary 
state in the sense of the topology defined by doo ■ Here doo stands for the infinity- Wasserstein 
distance on the space of probability measures (see section [s] for a definition) . We also show 
that if the repulsive-attractive potential W is singular enough at the origin, for example 
W{x) ~ —\x\^/b as |x| — )• with b < 3 — N , then the potential is so repulsive in the short 
range that solutions can not concentrate on an hypersurface, and this is independent of how 
attractive is the potential in the long range. To be more precise we show that for potentials 
with such a strong repulsive singularity at the origin, solutions can not converge with 
respect to the doo-topology toward singular steady states supported on hypersurfaces. 
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Whereas section |3] is devoted to instability results, section |4] is devoted to stability results. 
We show that if (C0)-(C2) hold with strict inequalities, then a radially symmetric solution 
of ([T]) which starts close enough to the spherical shell in the d^o topology will converge 
exponentially fast toward it. Under additional assumptions on the potential we can also prove 
convergence with respect to the da topology, a G [1, +00). In order for the stability results of 
section |4] to hold a certain amount of regularity on the solutions is necessary. Unfortunately 
weak L^'-solutions do not have this amount of regularity. This is why in section [5] we prove 
well posedness of classical C^-solutions. This covers a gap in the existing literature which 
mostly considers weak solutions. The results of section |4] are true for this class of classical 
C^-solutions. The aim of section [g] is to show examples of how to apply the general instability 
and stability theory in the case of power-law repulsive-attractive potentials: 

W{x) = ^ _ ^ 2-N <b<a. (3) 
a b 

For this family of potentials, conditions (C0)-(C2) can be explicitly formulated in terms of 
a and 6, therefore leading to an explicit bifurcation diagram for the stability of the spherical 
shell in . Finally in the last section, section [7| we perform numerical computations of 
radially symmetric solutions of ([T]) with power-law potential ([S]) and study their convergence 
toward spherical shell stationary state. Since a spherical shell is a highly singular function, it 
is challenging to perform such computations with traditional methods. This is why, following 
[MIOESESlEl], rather than simulating Q directly, we simulate the evolution of the inverse 
of the cumulative distribution of the radial measure associated to Since the inverse of the 
cumulative distribution of a spherical shell is a constant function, this approach has the virtue 
of smoothing the dynamics and this provides us with a robust numerical scheme. Our nu- 
merical simulations indicate the possible existence of integrable radial stationary states stable 
under radial perturbations in the parameter area corresponding to the fattening instability 
for power-law repulsive-attractive potentials, an issue that will be analysed elsewhere. This 
has already been proved in the particular case of 5 = 2 — and a > 2 in [34j . 



2. Preliminary section 
2.1. Radially symmetric formulation of the equation. 

Definition 2 (Radial Measures). We denote by V^{M.^) the space of radially symmetric 
probability measures. If jjL & V^i^^) then jl G "^([0, +00)) is defined by 

dfi{r) = / dii{x) and I djl{r) = / dfi{x) 

ri JTi<\x\<r2 Jo Jo<\x\<r2 

for all < ri < r2. We endow this space with the standard weak-* topology. 

Recall that 6ji G ^""(M^) stands for the spherical shell of radius R (see Definition [l]). 
The velocity field at point x generated by a spherical shell of radius R is given by vr{x) = 
—VW * dnix). Since W is radially symmetric, then by symmetry there exists a function 
uj{r, rf) such that 

vr{x) = -VW*5r{x)=uj{\x\,R)^^ (4) 

and one can easily check that this function oj is defined by see 181 for more details. Note 
also that if G P''(M^) then it can be written as a sum of spherical shells, ^ = Jq djl{rf), 



diu){r,rj) = -— I -^;-^{rei - r]y)da{y), (7) 
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and we conclude that 

I-+00 r+oo 

-{VW * fi){x) = - {VW * 6r,){x)djl{'n) = u{\x\,i])d^{ri) — 

Jo Jo pl 

Given T > 0, C([0, T]; P''(M^)) denotes the set of continuous curves of radial measures where 
continuity is with respect to the weak-* convergence. We say that fi E C{[0,T]]V^ (K^)) is 
a radially symmetric solution of ([T]) if /i S C{[0,T];V{[0, +oo)) satisfies the one dimensional 
conservation law: 

dtfi + drifiv) = (5) 

r+oo 

v{t,r) = / uj{r,ri)d(it{ri) , (6) 
Jo 

in the distributional sense. We will now give conditions for the velocity field to be well-defined 
by studying the properties of the function uo. 

2.2. Regularity of the function uj{r,ri). Let us remind that we assume that W is radially 
symmetric and belongs to C^(M^\{0}). The function uj{r, rf) defined by ^ is clearly away 
from the diagonal P = {(r, r) : r > 0}. Moreover, the derivatives of (jJ are given by 

1 r d^w 

(^N JdB{0,l) dxj 
and 

1 f fdW\ 
d2Uj{r, 7]) = — / V (rei - r]y) ■ y da{y), (8) 

(^N JdB(0,l) \oxiJ 

away from the diagonal. We need to investigate the behavior of u on the diagonal. Let us 
make the following definition: 

Definition 3 (Integrability on hypersurfaces) . A radially symmetric function g £ C(]R^\{0}) 
is said to be locally integrable on hypersurfaces if 

/ \g{x, 0)1 dx < +00 

J[0,l]^-i 

where x = {xi, . . . ,xi\f-i), or equivalently, if g{r)r^~'^ is integrable on (0,1) with g{x) = 
g{\x\). By an abuse of notation, we sometimes say g{r) is integrable on hypersurfaces. 

Lemma 1 (Regularity of the function oj). Let W{x) = /c(|x|) be a radially symmetric potential 
belonging to C3(M^\{0}). 

(i) If k'{r) is locally integrable on hypersurfaces then oj G C(M^). 

(ii) If k'{r), k"{r), and r~^k'{r) are locally integrable on hypersurfaces then uj £ C^(M^). 
(Hi) Suppose AW is negative in a neighborhood of the origin. Ifk'{r) is locally integrable 

on hypersurfaces but /SW = k" + {N — l)r~^k' is not, then for any R > 0, 

lim diuj(r,ri) = +oo. (9) 

Before proving the above lemma, let us discuss the result. Obviously the regularity of the 
function cj depends only on the behavior of W at the origin. Assume for simplicity that in the 
neighborhood of the origin, the potential is a powerlaw, that is W{x) = A;(|x|) = — 
for all X G -6(0, e), where b is possibly negative. Note that k'{r) < for r < e so the potential 
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is repulsive in the short range. Lemma [T] then claims that cj is continuous if 6 > 2 — and 
continuously differentiable if 6 > 3 — A^. Statement (iii) says that if2 — A^<6<3 — A^, then 
oj is continuous but its first derivative goes to +oo as (r, 77) approaches the diagonal. 

Finally, let us remark that(i) is sharp in the sense that the Newtonian potential is 
the critical one for the integrability on hypersurfaces. Precisely, Newton's Theorem asserts 
that the function oj associated to the Newtonian potential is discontinuous, it has a singularity, 
across the spherical shell. We now prove the Lemma: 

Proof. Let us prove (i). The function uj{r,ri) can be rewritten as 

u}{r,rj)=[ ei-VW{rei-y) ^r—rdaiy). 

JdB{0,r,) f^NV 

Seeing w as a function of x = rei and ij, we can apply Lemma [TO] from the appendix with 
:= dB{0,7]), <j)^{x) := {aNV^~^)~\ and G{x) := ei-VW{x). Since \G{x)\ is bounded by 
which is locally integrable on hypersurfaces, we obtain that to G C(M^). 
We now turn to the proof of (ii). It is simple to check that 



a a = k [r)^y^ +k{r)^-k (r)^ 
oxiOXj ^ '^•^ 



and then | g^,^ | is bounded by a radial function which is locally integrable on hypersurfaces 

given by a linear combination of k"{r) and r^^k'{r). Moreover, it has the regularity needed 
in Lemma 10, We now rewrite the derivatives diU!{r,r]) and (92w(r, r/) in Q and ([s]) as 

1 r d^W 
(7Nr]^~^ JdB{o,v) dxj 



diuj{r,ri) = - _ / -jr^{rei-y)d(T{y), 



1 r 



d2Uj{r, rj) = - — ^ / V ( ] (rei - y) ■ y da{y) . 



The reader can easily check that Lemma 10 applies similarly as before, so that w € C^(M^). 
Finally we prove (iii). Taking the divergence of (|4j) we obtain: 

(div vr){x) = -AW * 6r{x) = dM\x\ ,R) + {N- 1) "^^*^*;^^ (10) 

\x\ 

and therefore diuj{r,ri) can be written: 

diuj{r, 7]) = -{AW * <5„)(rei) - {N - 

r 

For < e < ro, let £ C°^(M+) be a cut-off function, such that = 1 on [0, e/2], and Xe = 
on [e, 00). Choose e such that the function —AW^{x) := —Xeix)AW{x) is nonnegative. Using 
Lemma 1 1 with r/i and 1^2 such that rji < R < r]2, and noting that dist (rei, dB (0,r/)) = |r — 
we find that 

lim -{AW * 6r^){rei) = +00. 



To conclude the proof, note that the functions (r, r/) 1— )• and {r,ri) 1— )• ([(1 — X£)A1^] * 

Sri){rei) are bounded in a neighborhood of {R, R). □ 
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2.3. The two radial instabilities. In this subsection we exhibit some elementary calcula- 
tions in order to understand under which conditions a spherical shell is a stable steady state. 
Rigorous results about stability and instability of spherical shell with respect to the transport 
distance will be provided in section 3 and 4. This subsection provide motivations for the 
rigorous results to come later. 

Definition 4 (Steady states). A probability measure /x G V{M.^) is said to be a steady state 
of the nonlocal interaction equation ([T]) if 

— {VW * = for all x E supp(//). 

We now show that if the attractive strength k'{r) of a repulsive-attractive potential W{x) = 
k{\x\) does not decay faster 1/r^ as r — )• oo, then there exists a spherical shell steady state. 

Lemma 2 (Existence of spherical shell steady states). Let W{x) = k{\x\) be a radially 
symmetric potential belonging to C^{M.^\{0}) and such that k'{r) is locally integrable on 
hyper surf aces. Let us assume that the potential is repulsive-attractive in the following sense: 
there exists Ra > such that 

k'{r) >Oforr> Ra, and k'{r) < for < r < Ra- 

Defining for r > 2Ra the function 

S(r) := inf k'{s) > 0, (11) 

r/2<s<2r 



we will further assume that 



lim r^Ti(r) = +oo. 



Then there exists at least a R > such that the spherical shell G T'(M^) is a steady state 
to 0. 

Proof. Note that from Q we directly obtain that a spherical shell 6ji £ V{M.^) is a steady 
state if and only if uj{R,R) = 0, that is, if and only if condition (CO) holds. Since k'{r) 
is locally integrable on hypersurface w G C(M^) due to Lemma [ij So the function F(r) := 
w(r, r) S C(M+). Using formula ([2]), we get 

f /" y — ^1 

F{r) = — / k'{r\y - ei\)- r • ei dcr{y) . 

(^N J 33(0,1) \y - ei\ 

Let us remark that (y — ei) • ei < for all y G dB{0, 1), and thus for 2r < Ra we easily 
get F(r) > 0. It is enough to show that there exists r > 2Ra such that F(r) < 0. In order 
to do this, we proceed as in [TT, Proposition 2.2] and divide the integral in the definition of 
F(r) into two sets: A := dB(0,l) n B{ei,Ra/r) and its complementary set A'^. Note that 
the integrand is positive on A and negative on A^. We will show that for r large enough the 
integral over the set A^ is greater in absolute value than the integral over the set A. It is easy 
to see that B := {y £ dB{0, 1) such that 2\y — ei| > 1} C ^4'^ and B ^ $ as soon as r > 2Ra. 
We first estimate the integral 



k'{r\y - eiDj^ ■ ei da{y) 

\y - eil 



\k'{r\y-e^\/^^—^da{y) 
|ei - y\ 



Let 6 be the angle between ei — y and ei and note that for ally £ A := dB{0, l)riB{ei, Ra/r) 
we have by the law of cosines 

ei-y a ^ 

■ ei = cos < 



\ei — y\ 2r 
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Using Lemma [s] from the Appendix with Ai = dB{0, 1) we then obtain 



k'{r\y - eiDj-!^ — ^ • ei da{y) 
|ei - y\ 



< ^ / \k'{r\y - ei|)|dcr(y) 
^ 5 T"^' l^'MI 1-^ n {|y - eil = s}\^N-2 ds 







7? rRa/r 
< / W{rs)\ s'^-^ds 

Jo 



2rN 70 ' ^ ^' - 

where we have used the fact that k'[r) is integrable on hypersurfaces to obtain the last 
inequahty. 

Since the integrand is negative in and since B C A"^ for r > 2Ra we have: 

k'{r\y-ei\)^ — ^■eida{y)< / k' {r\y - ei\) -f^ — -^-eidaiy). 
\y-ei\ Jb \y-ei\ 

Moreover, thanks to the law of cosines, j^Ef^ " = cos(— 0) < —1/4 for y & B, and then 
using (11 ) 

y — ^1 1 f 

k'{r\y - ei\)- r • ei d(j{y) < -7 / k' {r\y - ei|) da{y) 



<-\j^da{y)T.{r) < -C2 S(r). 



Condition ( |11[ ) on S(r) implies that C2S(r) > Ci/r^ for r large enough and therefore F{r) < 
for r large enough. Then the continuity of F implies the existence of a radius f > such 
that F{r) =0. □ 

The following proposition gives some hints about the stability properties of the spherical 
shell steady states. 

Proposition 1 (Instability modes by energy arguments). Assume that the radial interaction 
potential W is such that u £ C^{M.'^) and let be a steady state, that is uj{R, R) = 0. 

(i) If {CI) is not satisfied then by splitting the spherical shell into two spherical shells 
we can decrease the energy. More precisely there exists dr^ > such that, given 
< \dr\ < dro, 

E[{1 - e)6R + e6R+dr] < E[6r] 

if € is small enough. 

(a) If (C2) is not satisfied then by increasing or decreasing the radius of the spherical 
shell we can decrease the energy. More precisely there exists dro > such that 

E[6R+dr] < E[6r\ 

for all < |dr| < dr^. 
Proof. Let us introduce the notations 



E[l^, ^] := J // W{x - y)d^i{x)diy{y), 
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SO that E[fi,fj] = E[ij,], and 

11 r r 

W{rx — rjy) da{x)d(j{y) . 



E{r,rj) ■.= E[5r,Sr,] = ~ 



2ai 



N J JdB(0,l)xaB{0,l) 



Taking the derivative we get: 



dE 

dr 



r, r?) 



1 1 

2 7 JdB{0,l)xdB{0,l) 



1 1 



1 



"^CTN JdB{0,l) y^^N JdB{0,l) 
1 1 



VW{rx — rjy) ■ x da{x)da{y) 

VW{rx — rjy)da{y) ■ xda{x) 



2 CTn JdB{0,l) 



VW * 6r,{rx) ) • x da{x) 



-uj{r, 7]) . 



Since E{r,r]) = E(r],r), the Hessian matrix of E{r,ri) is given by 

H{r,rj) = - 



diuj{r,r]) d2Uj{r,r]) 
d2Uj{r],r) diU}{r],r) 



If 5ji is a steady state, i.e. uj{R, R) = 0, then VE{R, R) = and 

E{R + dr, R) = E{R, R) - ^diuj{R, R)dr^ + o{dr'^) 
E{R + dr,R + dr) = E{R,R) - ^{diuj{R, R) + d2Uj{R, R)) dr^ + o{dr^ 



(12) 
(13) 



The proof of (ii) follows directly from the Taylor expansion (13). By using the Taylor expan- 
sions (12) and (13) as well as the bilinearity of E[fj] = E[fi, fi]: 

E {1 - e)6QB(o^B) + ^^dB{0,R+dr) , (1 " (^)^dB{0,R) + ^^aB(0,R+dr) 

=(1 - €fE{R, R) + 2e(l - €)E{R + dr, R) + e'^E{R + dr,R + dr) 
=E{R, R) - -diuj{R, R)dr'^ - —d2Uj{R, R)) dr^ + o{dr^) 
from which (i) follows by taking e and drQ small enough. 



□ 



The following elementary Lemma shows that the instability condition diOj{R,R) > (i.e. 
(CI) is not satisfied) simply means that the divergence of the velocity field generated by 
the spherical shell is positive on the spherical shell. Being the velocity field "expanding", it 
makes sense that splitting the spherical shell into two reduces the energy as proven in previous 
Proposition [TJ 

Lemma 3 (Divergence of the velocity field). Assume the spherical shell 5r is a steady state, 
i.e., condition (CO); uj{R,R) = 0. Let vr he the velocity field generated by Sr, given by (Q. 
Then 



(div vr){x) = diuj{R, R) for all x G 95(0, R). 



(14) 



Proof. This is a direct consequence of (10) together with the fact that uj{R,R) = 0. 



□ 
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2.4. Well-posedness of L*^-solutions. Global existence and uniqueness of L^'-solutions of 
equation ([T]) was established in |1H Theorem 1] under some conditions on the interaction 
potential W: 

Theorem 1 (L^'-Well posedness theory). Consider 1 < q < oo and p its Holder conjugate. 
Suppose \7W G W^''^(M^) and fiQ £ LP(M^) n ^2(1^^) is nonnegative. Then there exists 
a time T* > and a nonnegative function fi e C([0, T*], LP(M^)) n C\[0,T*],W-^'P{R^)) 
such that ([1]) holds in the sense of distributions in X (0,r*) with /x(0) = fiQ. Moreover 
the second moment of x ^ IJ-itj^) remains bounded and the norm is conserved. Also, the 
function t — )• ||Ai(t)||f^p is differentiable and satisfies 

YAMt)\\%} = -iP-'^) f Kt,xrdivv{t,x)dx yte[0,T*]. (15) 
at j^N 

Furthermore, if ess sup AW < +00, then t — )• ||Ai(i)||^p does not grow faster than exponen- 
tially and we have global well-posedness. 

In the above theorem V2{R^) stands for the space of probability measure with finite second 
moment. We will refer to the solutions provided by the above theorem as L^-solutions. 

One can find in [1] that the authors extend the global-in-time well posedness L^-theory 
to repulsive-attractive potentials under suitable conditions. We summarize the result in the 
following theorem. 

Theorem 2 (Dealing with possibly growing at 00 attractive potentials). Assume that W{x) = 
A;(|x|) is a radially symmetric repulsive- attractive potential, W{x) = {Wr -\- Wji){x) = (A;_r -|- 
kA)i\x\) = k{\x\) with k G C^{{0,+oo)), Wa attractive (i.e. k'^ > 0), with VW £ W/j^(M^), 
1 < q < 00, and Wr compactly supported repulsive (k'j^ < 0). Furthermore, assume that k 
satisfies: 

(i) 35i > such that k"{r) is monotonic in (0,5i). 
(ii) 3(^2 > such that rk"[r) is monotonic in (0,^2). 
(Hi) D := sup^g(o,oo) WR{r)\ < 00 

(iv) There exists m such that is bounded and increasing. 

Then there exists a global in time solution for the equation ([T]) with compactly supported initial 
data Ho G LP(R.^), which is compactly supported for all t > 0. 

3. The Fattening instability and dimensionality of the steady state 

3.1. The radially symmetric case. This first subsection concerns radially symmetric solu- 
tions. We show that if the singularity of W at the origin is such that the kernel w is C^, and 
if condition (CI) is not satisfied, then a radially symmetric solution can not converge weakly 
toward the spherical shell stationary state. We also show that the same result holds if the 
singularity of W at the origin is so strong that the kernel u is not (and this is independent 
of how strong the attractive part of the potential is) . 

Theorem 3 (Instability of spherical shells: radially symmetric case). Let W{x) = k{\x\) be a 
radially symmetric potential belonging to C'^(M^\{0}) and such that k'{r) is locally integrable 
on hypersurface (so that 00 is continuous). Assume that the spherical shell 5r is a steady 
state, that is, (CO).' uj{R,R) = 0, and that one of the two following hypotheses hold: 

(i) k"{r) and r^^k'{r) are locally integrable on hypersurfaces (so that co is C^), and 
diUj{R,R) > 0. 
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(ii) AW is negative in a neighborhood of the origin and is not locally integrable on hyper- 
surfaces (in which case oj is not and lim {r,r])^v diU}{r,r]) = +00). 

Then it is not possible for an radially symmetric solution of ([5])-(|6| to converge weakly-* 
as measures to 6r as t —)• 00. 

To clarify the result, let us consider the case where the repulsive-attractive potential 
W{x) = A;(|x|) has its repulsive part described by a powerlaw. Let say, for example, that 
k{r) = -r^/h, for all r < 1 and k'{r) > for r > 2. If 2 - < 6 < 3 - iV, then AW is not 
locally integrable on hypersurfaces and therefore, according to (ii), whatever is the behavior 
of W{x) for |x| > 1, radially symmetric solutions can not converge toward the steady state. 
In other words if the repulsive singularity of the potential is equal to or stronger than \x\ 
then the potential is so repulsive in the short range that solution can not concentrate on a 
spherical shell, and this is independent of how attractive the potential is in the long range. 
On the other hand if 6 > 3 — then the kernel uj is In this case, the balance between the 
repulsive part and the attractive part of the potential dictates whether or not the spherical 
shell is an attractor: if diuo{R,R) > 0, then the repulsive part dominates and the spherical 
shell is not an attractor. 

We remind, see [591 [191 |l5] , that for 1 < p < cxd the distance dp between two measures z/, p 
is defined by 



(^p{^^P)= inf i / \x -y\Pd7r{x,y) 

where n(i/, p) is the set of those joint distribution functions with marginals ly and p. When 
p = +00 then the distance is defined as 

doo{y,p)= inf <^ sup \y-T{y)\ : T#p = v\. 



We now prove the Theorem. 

Proof. If conditions (ii) of the Theorem holds, then from ^ of Lemma[l| it is clear that there 
exists 6 > such that 

Vr/ G {R — 6, R + 6), r ^ u}{r,r]) is strictly increasing in {R — 6, R + 6). (16) 

Of course ( |16[ ) also trivially holds if condition (i) of the Theorem is satisfied. We proceed 
by contradiction. Assume that p{x,t) = p^t{x) is an L'p radially symmetric solution which 
converges weakly-* as measures to a spherical shell of radius as t — )• 00. 

Step 1. Assume first that pt converges toward 5r not only weakly-* as measures but also 
with respect to the doo-topology. This implies that the support of the radial solution pt to ([S]) 
converge to the point {R}. Choose T > such that supp(/it) C {R — 5,R + 5) for all t > T. 



Using the monotonicity property (16) we obtain that for t >T and for R—6 < ri < r2 < R+6 

j-R+S 

v{t,r2) - v{t,ri) = / uj{r2,r]) - uj{ri,ri) dpt{ri) > 
J R-5 

where i) is the velocity field in radial coordinate defined by ([6]). Therefore for all t > T the 
function r — )• v{t,r) is increasing on {R — 5, R -\- 6). Let ri{t) and r2{t) be two solutions of 
the ODE r'iit) = v{t,ri{t)), i = 1,2. Since 

^(r2(t) - n{t)f = 2(r2(t) - r,{t)){v{t,r2it)) - ^)(t,ri(t))) > , 
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we easily see that if for some time t > T, ri{t) and r2{t) are in (R — 6, R + 6), then their 
distance increases. This contradicts the fact that the support of fit is converging to the 
point {R} as t — >• oo. Let us be more precise. Since fix is supported in {R — 6, R + 6) 
and is absolutely continuous with respect to the Lebesgue measure, there exists i?i and 
R2 in {R — 6, R + 6), Ri 7^ R2, such that fiT{x)dx = 1/3 and fiT{x)dx = 1/3. 
Consider the ODEs r[{t) = v{t,ri{t)), ri{T) = Ri, i = 1,2. Clearly ri{t) and r2{t) remain in 
{R — 5,R + 5) for alH > T (otherwise the support of would not stay m. {R — 5, R + 5)) . So 
|r2(t) — ri{t)\ > \R2 — Ri\ for all t > T and the support of fit can not converge to the point 
{R}, which contradicts our assumption. 

Step 2. Assume now that [it converges weakly toward 5r but does not converge with 
respect to the doo-topology. From the continuity of the function r] uj{r, rf) together with ^ 
it is clear that v{r,t) converges pointwise to uj{r,R) as t — )• 00. Since the support of /tj does 
not converge to the set {R} there is a sequence of times at which there is always non-zero 



mass in (0, i? — e) U + e, +00). Since ijj{R, R) = 0, the monotonicity condition (16) implies 
that (jj{r,R) < for all r S {R — €,R) and Lo{r,R) > for all r G {R,R + e) as long as 
e < 5. Because of the pointwise convergence of v there exists a time T > such that for all 
t > T, v{t, R — e) < and v{t, R + e) > 0. So after this time T mass cannot enter the region 
[R — e, R + e]. This together with the existence of a time t > T for which there is some mass 
in the complementary of [R — e, R + e] contradict the weak convergence towards 6r. □ 

Remark 1. In Step 2 of this proof, since we are dealing with radially symmetric solutions, 
the problem is essentially one dimensional and the characteristics are ordered. This allows 
us to exclude the possibility of an U' solution converging toward a spherical shell even if this 
convergence is very weak and the support of the solution does not converge. In the non radially 
symmetric case we will be only able to exclude convergence in d^o . 

3.2. The non-radially symmetric case. In this subsection we consider non-radially sym- 
metric solutions and we investigate whether it is possible for an LP-solution to converge toward 
a steady state supported on an hypersurface which not necessarily a sphere. Indeed in numer- 
ical simulations [SOlEQ], it is observed that depending on the choice of the repulsive-attractive 
potential W, solutions of ([T]) can either converge to steady states which are smooth densities 
or to singular steady states which are measures supported on an hypersurface. We consider 
steady states p, of the form 

fix)dfi{x) = [ f{x)<P{x)d<j{x) V/ G C(M^) (17) 
Jm 

where is a compact hypersurface and da is the volume element on Ai. Roughly 
speaking, we prove that if the the potential is as singular or more singular than at 
the origin, then it is not possible for an L^-solution to converge toward such a steady state 
with respect to the doo-topology. We also prove that the same result holds if the potential is 
less singular than and if the divergence of the velocity field generated by such steady 

state is strictly positive on its support. 

Theorem 4 (Instability of Spherical Shells: Nonradial case). Let W{x) = k{\x\) be a radially 
symmetric potential which belongs to C^(M^\{0}). Assume that \\m.r^Q KW {r) = —00 and 



that close to the origin AVF(r) is monotone. Let jl be a steady state of the form (17) with 
M. being a compact hypersurface and let v be the velocity field generated by Jl, that is 
V = —VW * pL. If one of the two condition holds: 
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(i) AW is locally integrable on hyper surf aces, (j) £ L°°{M) and 

(div v){x) := —{AW * ll){x) > for all x G supp /2, (18) 

(ii) AW is not locally integrable on hypersurfaces and (/>(x) > (/>o > for all x £ Ai, 

then it is not possible for an solution of ([T]) to converge to jl with respect to the doo-topology 
as t ^ CO. 

Before to prove this Theorem, let us make some remarks: 

Remark 2. According to Lemma^ the result of Theorem^ (i) of the previous subsection 
can be reformulated as follows: assume that the spherical shell 6r is a steady state and let vr 
be its velocity field. If 

(div vr){x) = -{AW * 5r){x) > for all x G 55(0, R) 

then it is not possible for an LP radially symmetric solution to converge weakly-* as measures 
to Sr as t ^ oo. So conditions (i) of Theorems^ and [4] are essentially the same. Similarly 
condition (ii) of both Theorems are also essentially the same. In this sense Theorem^ can be 
seen as a generalization of Theorem^ to the non-radially symmetric case. 

Remark 3. The assumption liuir^o AW{r) = — oo simply guarantees that the potential W is 
strongly repulsive at the origin. The monotonicity of AW{r) in a neighborhood of the origin 
is not essential to the proof and could be replaced by weaker hypotheses. But in practice all 
potentials of interest satisfy this monotonicity condition. 

Remark 4. Part (ii) of the Theorem, roughly speaking, states that if the repulsive- attractive 
potential W is more singular than at the origin, then whatever is its attractive part, 

it is not possible for an solution of ([T]) to converge with respect to the doo-topology toward 
a singular steady state supported on an hypersurface. So we see that the dimensionality of 
stable steady states depends on the degree of singularity of the potential. For such potential 
with a strong repulsive singularity at the origin, steady states are expected to be absolutely 
continuous with respect to the Lebesgue measure. 

Theorem |4] is a direct consequence of the three Lemmas to follow. 

Lemma 4 (Approximating the divergence of the velocity field). Let W be as stated in The- 
orem |4] and let p, be a compactly supported probability measure not belonging to L^. Suppose 
there exists a Holder continuous function AW > AW such that 

- AW * fl> on supp(/i) , (19) 

then it is not possible for an solution of ([T]) to converge to p, with respect to the d^o-topology 
as t ^ oo. 

Proof. We proceed by contradiction. Let fit be an LP solution such that liuit^oo rfoo(/Ut, /i) = 0. 
We are going to show that there exists a T > and an e > such that for alH > T 

{AW * fJ.t){x) < -e for all x £ supp(^t), (20) 



and combined with the equality (15): 

d 



^^„/^t|lLP = (P- 1) / _ {AW * fit)lJ'tdx, 
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this guarantees that a subsequence fit„ converges weakly in to an function, which 
contradicts the assumption that p, ^ L^. Let us prove (20). Write 

AW *nt = Ary * {nt -fl) + AW + {AW - AW) * (21) 

and note that since AW > AW the third term is negative for all t and x. Since AW is 
continuous so is AW * fl and, therefore, ( 19 ) implies that there exists an e > and an open 
set Q containing the support of p, such that AW * p < on Q. Here we used that the 
supp(/Z) is a compact manifold. Note that since limt_>oo doo(/^ti A) = the support of fj^t will 
eventually be in To estimate the first term of (21), we consider Tt '■ - 
pushing forward /it to n, i.e. Tt#^it = A- Then 



pN 



a map 



AW * {fit - Jl){x) 



< 



\AW{x - y) - AW{x - Tt{y))Wt{y) 

c\y -Tt{y))fdfit{y), 
since this inequality is true for any map % pushing forward fit to 



< 



\\Aw*{f,t-mu 

so that for t >T with T large enough, 

\\AW * {fit - 



< C doo(/Ut,/i)^, 



< 



Since AW > AW, the last term of (21) is negative, so that AW * fit < 0. 



□ 



In order to conclude the proof we now need to show that under the hypotheses of the 



theorem there exists a Holder continuous function AW > AW satisfying (19). Define 

AW'{x) :-- 



AW{x) if \x\ > e 
AW{eei) if \x\ < e 



The function AW^{x) is obviously Holder continuous and, due to the monotonicity of AW 
around the origin we have AW' > AW for e small enough. We are left to show that AW'*fL < 
on the support of fl for e small enough and this is done in the following two Lemmas. 

Lemma 5 (Continuity in e of the divergence of the velocity field). Let fl, Ai and W be as 
stated in Theorem[4] (i). Then AW"*!! converges uniformly on A4 toward AW* p. Therefore 
there exists e > such that AW^ * ft < on the support of p,. 

Proof. Since A4 is and compact, using Lemma [s] from the appendix, there exist constants 
6,Ci,C2> so that 



Ci / g{r)r^-'^dr < 



MnB{x,e) 



g{\x - y\)da{y) < C2 I g{r)r^~^dr 



for all X G A^, for all e < 5 and for all nonnegative function g locally integrable on hyper- 
surfaces. Since AW{x) is radial and goes to — oo monotonically as \x\ — )• 0"*", we clearly have 
that AW{eei) — AW{x — y) > for all y E B{x, e) if e is small enough. Then we obtain that. 
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for all X e M, 



\{AW'*fi){x) - {AW*fl){x)\ 



< Co 



/MnB{x,€) 
< C2 \\4>\\l°-{M) 

and we conclude using the fact that AW is integrable on hypersurfaces. 



AW{eei)-AWix-y) (t>{y)dcj{y) 



r'^-^dr 



□ 



Lemma 6. Let ^, M. and W he as stated in Theorem [4] (ii). Then there exists e > such 
that AW^ * fL < on the support of p,. 

Proof. Choose rg as in Lemma [s] and also small enough so that AW{r) < for all r < rg. 
For e < ro we then have 

{AW'*fi)ix)< f AW{x-y)dJi{y)+ I AW{x - y)dji{y). (22) 

Je<\x-y\<ro J\x-y\>r(, 

Since AW is bounded on 5(0, 2 diam(A^)) \ i?(0,ro), the second term is uniformly bounded 
for X G . We use Lemma [S] to estimate the first term of ( 22 ) : 



'^<\^-y\<ro 



AW{x - y)dfi{y) 



Mr\{y:t<\x-y\<ro} 



AW{\x-y\)(l){y)da{y) 



< 



ro 



AW{r) \M n dB{x, r)\^N-2 dr < C^o I AW{r)r^-^dr^ 



and since AW is not locally integrable on hypersurface and AW < on [0,ro], the last 
integral goes to —00 as e — t- 0. Then, for e > small enough, AW'' * /i < on supp {fi). □ 

4. Stability for Radial Perturbations 

In this section, we give sufficient conditions for the stability under radial perturbations 
of Sji stationary solutions in transport distances for the system ([5])-(|6]). Let us denote by 
7^2(1^'^) the set of radial probability measures with bounded second moment. 

Here, we will work with radial solutions with the following hypotheses of minimal regularity 
(HMR): we assume that for any given fiQ £ P^(M^), there exists n £ ^C([0, T], P^'(I^^))> 
with fj,t = fiQ for t = 0, such that 

»+oo 



v{t, r) 



uir,r])diltir]) € L'{{0,T) xR^)) 



for all T > and their corresponding radial measures fit satisfy ^ in the weak distributional 
sense. Moreover, they satisfy that f^rdjlt{r) is an absolutely continuous function in time 
for which 

rdfit{r)= / v{t,r)djit{r) (23) 



dt 

holds a.e. t > 0. Furthermore, if fiQ is compactly supported, we assume that 

ri{t) = minjsupp (/it)} and r2(t) = maxjsupp {jit)} , 
are absolutely continuous functions with ^ri{t) = v{t,ri{t)) a.e. t > 0, i = 1,2. 
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The existence theory developed in Section [5] ensures that smooth classical solutions satis- 
fying (HMR) exist for /ig £ P2(^^) ^^'°°(^^) initial data under suitable assumptions on 
the potential. Therefore, we assume in this section that our radial solutions satisfy ([5])-(|6]) 
with oj given by ^ verifying suitable hypotheses specified in each result. 

Theorem 5 (Stability for local perturbations). Assume co £ C"'^(M^) as given in ^ and 
that 6ji is a stationary solution to (jsjl-Q, that is, the condition (CO); uj{R,R) = 0. Let us 
assume that (CI) and (C2) are satisfied with strict inequality, that is: 

diU}{R, R) <0 and diu;{R, R) + d20o{R, R) < . (24) 

Then there exists > such that if the initial data G (I^^) satisfies supp(/lo) C 
[R — Eq, R + Eq], and for any solution to ([5|-([6]) with initial data satisfying (HMR) we get 

d2{k.5R)<Ce-^\ 

for any < 7 < — max {diOj{R, R), j^u}{R, i?)) for suitable C. 

Proof of the Theorem. Since we have assumed that the solutions to ([5])-(|6]) satisfy the regu- 
larity conditions (HMR), then T{t) := diam (supp(/it)) = r2{t) — ri{t), and 



e{t) := / rdfit{r)-R, 
Jo 

are absolutely continuous function of t > 0. We will proceed by contradiction. 

We define T := min{t > 0; T{t) + \Q{t)\ > 4eo} and let us assume that T < 00 for all eo > 
close to 0. Note that T > by continuity of T{t) + |0(t)|, since supp (/xq) C [R — Eq, R + Eq] 
implies that 

i-R+eo 

r(0) + |e(0)| < 2eo + / \r- R\ d/io(r) < 3eo • 

JR~£0 

Now, for t G [0, T], supp(/it) C [R- Aeq, R + 4eo], since 

Init) -R\< \r,{t) -{R + Q{t))\ + |e(t)| < (r2(t) - ri{t)) + \Q{t)\ 

= r(t) + |e(t)| <4eo, (25) 

using that the center of mass Q{t) + R is obviously in [ri{t),r2{t)], for all t > and the 
definition of T. 

Then, for t £ [0, T], Taylor expanding to order one and using that dioj is uniformly contin- 
uous on [R-4:Eo,R + 4eo]^ together with (HMR), we get 

d d d f°° 

jTit) = -r2{t) - -ri{t) = v{t, r2(t)) - v{t, n{t)) = / [a;(r2(t), 7?) - t^(ri(t), r/)] dfit{r]) 



[dMn{t),r,){r2{t) - ri{t)) + g{ri{t),r2{t),7^)] dfniv), 
where g satisfies 



hm ( sup ^(lllI^Uo. (26) 

\r2-ri\^0 \ne[rur2] 1^2 -^1 / 



Since (26) is satisfied, the integral of g can be estimated as follows 

g{ri{t),r2{t),ri)dfit{ri) = / g{ri{t),r2{t),r]) dfit{r^) = o{r2{t) - ri{t)) = o{V{t)). 
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Proceeding with the same argument as before using (26), we can estimate 
-r{t) = {r2{t) - n{t)) dMn{t),r,)d(it{ri)+o{T) 

[diuj{R,R) + {dMn{t),n) - dMR,R))] dfit{ii) + o{T). 



Since ry G supp(/Uf) = r2(t)] <Z [R — A£q,R + Aeq] thanks to (25), we can then use the 

uniform continuity of diuj on [R — 4eo, R + 4eo]^ to get: 

\diuj{n{t),r]) - diOjiR, R)\ < C\ri{t) - R\ + \v - R\ < C\n{t) - R\ + \r2it) - R\, 



for any r/ G supp(/it). We can then use (25) again giving 
d 



dt 



T{t) = OMR, R){r2{t) - ri{t)) + o{T) + o(|e|). 



On the other hand, we can also estimate using (23) 
d 



dt 



v{r,t) djlt{r) 



JO 



a;(r,r?) djj.t{r) dfitif]) 



JQ 



[uj{T],r]) + diu}{7], r]){r - rj)] dfit{r) dfit{r]) + o{T) + o{\Q\), 



where we have again used an argument as in (26) to estimate the rest term of the Taylor 
expansion, and we use it once again to obtain 



d_ 
dt 



oj{R, R) + -^u;{R, R)i7] - R) + 9it^(r?, r/)(r - r/) 
Jo L "-n- 

+ o(r) + o(|e|) 

^oo{R, R) (^j^ r/ dfitiv) -R^+ dMR, R) (^j^ r dfi ' 

+ o(r) + o(|e|) 
^..(i?,i?))e + o(r) + o(|e|). 

We now combine the estimates on F and to get: 



dlit{r) dfitiri) 



d 
dt 



(r + |e|)(0 <max(5iw(i?,i?), ■^uj{R,R) ] {T + \e\) (t) + o {F + \e\) . (27) 



Let us point out that all the o(T + |0|)-terms can be made uniformly small in the interval 
[0,T] by taking eq small by their definitions and using that supp(/ij) C [R — 4eo,-R + 4eo] 
in [0, T]. More precisely, let 7 E (0, — max (9ia;(i?, i?), jj^uj{R, R))). We can choose eo > 
small enough for the rest terms of (27) to satisfy: 

o{m + \em 



m + \m\ 



< 



max ( diUj{R,R), -^lj{R,R) 
dR 



7> 



(28) 



for any T{t), Q{t) since r(t) + |e(t)| < Aeq for ah t G [0, T] due to i^. Then ([28]) is satisfied 
for all t £ [0, T], and thus, 

d 



dt 



(r + |e|)(t)<-7(r + |e|)(i), 
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SO that for t G [0,r], 

(r + |G|)(t) < (r + |e|)(o)e-^*. (29) 

In particular, for any time t £ [0,T], (r + |e|)(t) < (r + |e|)(0) < 3eq and thus, using 
the continuity of (r + |0|)(t) since T < +00 we can continue up to T > T satisfying 
(r + |G|) (t) < 4eo contradicting the definition of T. Thus, T = 00 for small enough eo 



and (29) then holds for all t > 0. Thanks to (25), this implies the exponential convergence of 



d2{fi'tjR) to 0: 

/•oo 

d2{i^t, Snf = / {r-RfdfLt{r) < max {\n{t) - R\'', \r2it) - R\^) < (T + \e\f (t) < Seoe"^* 
Jo 

for all t > 0. □ 

Remark 5. Lemma [l] gives sufficient conditions to get the assume regularity uj G C^(M^). 
Previous Theorem holds for all radially symmetric potentials W{x) = k{\x\) belonging to 
C^(M^\{0}) such that k"{r) and r^^k'{r) are integrable on hypersurfaces. This applies also 
to the next result for nan local perturbations. 



Remark 6. The first part of condition (24) implies intuitively that the velocity field created 
by given by uj{r,R) is decreasing at r = R and therefore, particles are pushed locally in 
space and in time towards radius R for small perturbations. 

From now on, we denote by f{t, •) the pseudo-inverse of the distribution function of the 
radial measure fit, that is 



(/p(t,0 = mf|rGM+; d/it > ^ j • (30) 

if then satisfies 

/■oo 

dMt, e) = v{t, ^{t, 0) = / w(v9(t, 0, r/) dfitiv)- (31) 
Note that by the definition of ip, 

[ dfitiv) = [ dC. (32) 

J[ri,r2] J {S.;r-i<ip{t,S,)<r2} 

In the next theorem, we will work with solutions to system ([5])-([6]) satisfying (HMR) for 
which the pseudo-inverse of the distribution function is an absolutely continuous function on 



time satisfying (31) in the classical sense a.e. in t. Solutions obtained in Section [5] do satisfy 
these conditions. 

Theorem 6 (Stability: Tail control). Assume oj £ C^(M^) and that 6r is a locally-stable 
stationary solution to that is, ui{R,R) = and the local stability condition (24) 

holds. Assume moreover that the velocity field associated to Sji verifies 

uj{r, R)>0 on (0, R), u}{r, R) < on {R, 00), and diLo{0, R) > , 

and the following long-range controls on the interaction potential: for some a > 1, there exists 
A > such that 

oj{r, v)< \ - Ar" for (r, r]) e R+ x [R - X, R + X], (33) 
A 

sup\dM;ri)\<\a + Vn, (34) 
[o,A] A 
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|w(r,?7)| < -(l + r")(l + ?7°) for {r,rj) G M^. 
A 



(35) 



Then, for any solution to satisfying (HMR) and ( |31[ ) with initial data /io S P2(^^) 

such that /io({0}) = 0, and da{fiQ-,6ii) is small enough, 

lim da{fit,SR) = 0. 

I— >oo 

Remark 7. // we assume that the initial condition is compactly supported, then the long- 



range controls (33), (34), (35) on the interaction potential are not required anymore. Those 



are only necessary to control the behavior of the tail of the distribution and its interaction 
with the rest. 

Proof of the Theorem. 

Step 1.- "Claim: Given fi G 7^2(^^)- If daif^j^n) is small, then the associated velocity 
fields to fi and 5r share some confining properties": For any > small enough, thanks to 
our assumptions on uj, we can show that there exists A > such that if da{fi, 6ji) < A, then 



v{r) > Ci r > on (0, 
v{r) > vi > on ['d,R--d], 
v{r) < —vi on [R + 'd,oo), 



(36) 



where v{r) is the velocity field associated to /i by ([6]). To prove the first inequality, notice 
that da{jl,8ii) < A implies that 



dfi{ri) < K-'^''^ |r/-i?|"d/i(r?) < At (37) 

is small. We can then estimate the velocity field v for < r < < A: 



v{r) 



U}{r,r])dfi{ri) = / [uj{0,ri) + rdiU}{e,r])] dfi{r]) 
Jo 

r / diuj{9,ri) dfi{r]) + r / diuj{9,ri)d/l{r]). 

J[R-yfK,R+y/K] J[R-^^K,R.+^^XY 



Note that a;(0, rf) is equal to zero by definition. We then use (34) to get the following estimate: 



[R-y/k,R+y/KY 



dioj{9,r])dfi{r]) 



J[R-y^,R+y^Y 

<c(Af +d„(/i,5fi)-) <CA"/2. 



Now, if A is small enough and r € [0,??], then thanks to (37) and an argument as in (27) we 
conclude 

v{r) > rdiu{0, R)(l - [ dfi{r]) - CrA^^A - CrA"/^ > 

\ J[R-^,R+^Y J ^ 
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The second inequality (36) comes directly from assumption (35) and the continuity of u: 



for r £ [djR — and A small enough, 



v{r) 



[R-^/A,R+^/A] 



> (oj{r, R) - \/A||(92w| 



uj{r, rf) dfj.{ri) + 



uj{r,ri)dfi{r]) 



C[ (l + r/°)(iA(r/) 

J[R-^/A,R+^/A]'' 

>^a;(r,ii)-CA"/2, 



[R-^/A,R+^/A] 



where we have used (37). Since uj{-,R) > on {0,R), if A > is small enough, v{r) > on 
['d,R-^]. 

For the last inequality in (|36|), we can write the velocity field as 



[R-y/A,R+^/A] 



u{r,r])dfL{r]) + 



[i?-VA,i?+VA] = 



uj{r, 7]) dfi{r]) 



(38) 



and estimate the second term of (|38|) using (|35|) and (p37|) to obtain 

(1 + r/") (i/l(r/) < (:7A"/2(i + 



[R~^/^,R+^/^Y 



uj{r, ■n)dfl{ri) 



< C7(l + r") / 



[R-y/A,R+y/AY 



Let us distinguish two cases. In the set ^ (x ~ l)) ° which is equivalent to ^— Ar" < —1, 
we deduce that there exists Ci > such that ^ — Ar" < — Ci(l + r°). We can then control 
the first term of ( 38 ) using ( 33 ) and ([37]) to get 



L 



[iJ-VA,_R+VA] 



a;(r, r/) djl{r]) < ( ^ — Ar 



[/?-Va,/?+Va] 



dfi 



<-Ci(l + r-)(l - / djA <-^(l + r-) 

\ J[R-^,R+y^Y J ^ 



For r £ I : = 



.A VA 



1)) 



1/a 



, we use the assumption that u}{-,R) < on the compact 
interval /. By continuity of w, we thus have that for A > small enough and r £ I, 

max|a;(r,?7); r £ I,rj £ [R- VA, R + VA]^ := -C2 < 0, 

and thus. 



uj{r,r]) dji{vi) < -C2 



dpi<~< -C73(l + r") <0, 
[R~^/A,R+^/A] ^ 



'[R-v^,i?+VA] 

for r £l and A small enough using ( |37[ ). Then, ( |38| ) becomes 

v{r) < min(Ci, C3) + CVA^ (1 + r°) < -C4(l + f 
for any r>i? + ??, ifA>Ois small enough. 
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Step 2.- "Claim: If fiQ is close enough to Sr, then jit satisfies (36) at all times." Let ip{t,^) 
the associated pseudo-inverse function associated to jit by ( 30 ) . We assume that /xq satisfies 

For any e > 0, we can estimate |(/?(0, \/e) — R\ as follows: 

£: > dai^R, jio) > di{6R, jlo) = f \lp{0,0 - R\d^ 

Jo 

> max 1(^(0, - R\ d^, ^ \ip{0, - R\ d^^ . 

Since if is not decreasing, if ip{0, ^/£) < R, then \ip{0, S.) — R\ > \ip{0, -y/e) — R\ for ^ G [0, ^/e\. 
If 99(0, Ve) > R, then 193(0, - ^| > Iv^lO, y/e) - R\ for ^ E 1], so that 

e > min (Vi|(/p(0, - (1 - \/i)|(^(0, ^i) - , 

which provides the estimate \(p{0, y/e) — R\ < ^/e. Similarly, |(/9(0, 1 — y/e) — -R| < y/e, so that 

ip{Q,[^/e,l-^e\)(Z[R-.fe,R + ^/e\. (39) 

Let us define Te{t) := (^(t, 1 — y/e) — ip{t, ^/e) and 



Notice that for ^ G 1 — -y/e] , 

\ip{t,0-R\<re{t) + ee(.t). (40) 

For < e < eo < A, we define T := min{t G [0,r]; re(t) + \@e{t)\ > £o or da{5R,jLt) > Eq}- 
Thanks to ( |39[ ), T > by continuity for < e < min(A, sq) small enough. We will show that 
there exists e > such that T = +oo. Assume by contradiction that T < oo. By definition 
of T, we have dai^R, jit) < eo < A for t G [0, T]. 

Thus, jit satisfies (36) for i G [0, T], v{t, •) is positive on [0, i?— ??] and negative on [R+'d, oo), 
and then (31) implies that ^ G [0, 1], \ f{t,S,) — -R| < max(|(/?(0, ^) — -R|,i?). In particular, by 
(37) we get 

/ \ipit, e) - i?r < / mo, o - Rr + n dc < c^e. m 

For t G [0,T], we deduce that 

jTsit) =v{t, ifit, 1-V^))- v{t, ifit, V^)) 

= / [^(93(i, 1 - \/e), ??) - \/e), djLt{ri) 

+ / ['^{'P{tA-V£),v)-^{'P{t,V^),v)] djitiv)- 
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The first term can be estimated as it has been done for T{t) in the proof of Theorem [5] To 
estimate the second term, we use (35), (40), and (41) to conclude that 



/ 



[uj{ip{t, 1- ^/e),r|)- uj{ip{t, ^/e), ??)] dfit{r]) 



< ^ [1 + min {ipit, v^r, ^{t, 1 - v^r)] [ 



[<p(t,V^),95(t,l-v^)]- 



< ^[i + o(r, + |ee 



CV^+ [ \^it,C)-R\^dC 

J[v/i,l-v^]- 



<cVi + o(re + |e,|). 

The same can be done for 0^, and we obtain 

^(r, + |e,|)(t) <max(^aia;(i?,i?), ^oj{R,R)^ (T, + |e,|) (t) + o (F, + |e,|) + C7^. 

As it has been done in the proof of Theorem [5} eo can be chosen small enough such that this 
implies for t G [0, T], that 

I (r, + |e,|) it) < -7 (r, + |e,|) (t) + cVe, 

where 7 := 5 \max {diUj{R, R), ^uj{R,R))\. Then, for t € [0,r], 

(r, + |e,|) (t) < max (^(r, + |e,|) (o), ^Vi) < c^/i, 



and 



da{SR,fltr<[irs + \Qe\){t)r+ [ \^{t,0-RrdC 

<max(^(re + |ee|)(0),^Vi) + C^i < C^i, 



due to l\3m and (41). 



If e > is smah enough, this implies that for t £ [0,T], {Ts + \@£\){t) < ^, and 



dai^R, IJ-{t)) < ^. By a contradiction argument similar to the one used in the proof of 
Theorem [sj this shows that if e > is small enough, then T = +c«, and ( |36| ) is satisfied at 
all times. 



Step 3.- "Claim: Asymptotic convergence of jit to 5r:" Since (36) is satisfied for all t > 0, 
v{t, r)>Cir on [0, d], and then ip{t, C) > v?(0, ^e'^^* if (fit, < i9, due to We can thus 
estimate, using (32): 



d^ 



f{l + n dfit{r) <(1 + r dMr) = (1 + [ 

Jo Jo J{i;if 

<(l + i?") / dC = {l + Rn djloir). (42) 

J|f;ip(0.f)<i9e-<^in Jo 



Since (36) is satisfied, we claim that v{t,r) > vi on ['&,R — '&], and then tf{t,£,) G [0,R — •&] 
for t implies that (^(t - ^,0 < ^- 
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To see this, we make use of (|36]) to get v{t,-) > on [0, iJ-i?]. If (fit, £[0,R-^], (/?(•, 
is thus increasing on [0, t]. If (p{t- ^,C)> then ip{[t - ^,t],(,) C [&, R - i?] 



R 

R><f{t,C)-^{t--,0= / 
^'i Jt 



n R 

Vl 



which is absurd, thus (p{t — ;^,^) < "(? as desired. We can then estimate, for i > ^, using 



(32), 



■d 



R-i) r-R-d / r-d r-R-d 

(1 + r") dfitir) < (1 + ii") / dfit[r) = {1 + R^) ( dfitiO + / dfitiO 
JO \J0 

V-^{?;'/'(t,?)6[o,i?]} J-& J 



< + / + + / 

= (1 + i?") / (iAo(r) + (1 + Rn / (i/it-i?M(r). 

JO Jo 



dC 



Now, using a similar argument as in (42) since (f{t — |^,^) < in the last integral, we get 

{l + r")dfLt{r)<2{l + R") d^r). (43) 

Jo 



Thanks to ( 36 ) , if ^) > + i9, then 



= <^(0, + / Hs, ^{s, 0) ds < ^{0, - vit . 
Jo 



In particular, (^(t, < (/?(0,^)" and, thanks to (31), we get 

{l+r^)d(it{T)= [ il + ip{t,Cr)dC< [ 

J U 



R+'ff 



{1+ ifio, ^r)dc 



R+i)+vit 



(l + r-)d/io(r). 



Let e > 0. Thanks to (43), (44), there exists r > such that for any t>T, 

(l + r")dAt(0 < V^- 



R-t 



(44) 



(45) 



Then, in particular, for any t > t, djit{r) < \fe and dfLtir) < that is 

ip{t, [y/£,l- y/e]) C [R-'d,R + ^] and (Te + \Qe\) (t) < 3i? due to (go]). Now, for t > r, 
with an argument similar to the one used in Step 2, we get 



dt 



ir, + \Q,\){t)<max(dMR,R), (r. + |e,|)(t) + o(r, + |e,|) + C7V^. 

(46) 
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Since for any t > t, ip{t, takes its values in the compact set [i? — i? + inde- 

pendent of e, we can apply an argument similar to the one used in Step 2. Choose "i? > such 



that (46) implies 



I (r, + |e,|) (t) < -7 (r, + |e,|) (t) + cv^, 

and then, there exists some T >t such that for t >T, 



(r, + |e,|)(t)<2^. (47) 

7 



To conclude, we notice that thanks to (45) and (32), 

/ \ipit,o-Rrdc< I max(^Mv9(t,e)-i?r)^ie 



< 2i9"Ve+ / \r - Rl'^dfitir) 



which, together with (47), implies that for t >T, 

[ \^^t,^)-Rrdc+ [ \ip{t,c)-Rrd^ 



da{h^^R) 



Since this is true for any e > 0, it shows that daifj-ti ^r) — ?• as t — )• oo. □ 

5. Existence theory 

Existence and uniqueness of weak solutions for the aggregation equation in 7^2(1^^) H 
LP(M^) have been given in [4Ul 13 [TT]. Weak measure solutions to the the Cauchy problem 
for the aggregation equation ([T]) where given in [18] under the condition that the potential is 
smooth except possibly at the origin, the growth at infinity is no worse than quadratic, and 
the singularity at the origin of the derivative of the potential is not worse that Lipschitz. This 
section is aimed to give an existence theory of classical solutions for the aggregation equation. 

We will denote by ^"•^(M^), 1 < p < oo and m G N, the Sobolev spaces. 

Theorem 7 (Existence of classical solutions). Let W satisfy 

VW G L\R^), D^W G lI(M^), iAW)+ G L~(M^) (48) 

Then, for any initial data po{x) G W^'°°(M^), there exist classical solutions p G C^([0,T] x 
R^) n WI'^{R+, W^'~(M^)) to (0. Moreover, if po{x) G W^'°°(M^) /or k G N, k > 2, then 
p G C^((M^) X [0,r]) r\Wi~^'°°{R+,W~^{R^)). Furthermore, assuming in addition that 
Po G L^{'K'^) with bounded second moment, the solution is unique. 



Proof of the theorem. Step 1: A priori estimates. In this step we assume that the solution is 
smooth as needed. This assumption will be removed in the next step. 
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We consider first x G such that p{x,t) = \\p{t)\\^. Then, Vp{x,t) = 0, and 



dtp{x,t) = Vp{x,t){VW * p){x,t) + p{x,t){AW * p){x,t) 
<p{x,t){{AW*p)+{x,t)) 
< ||(AW^)+|Loo p{x,t), 

so that 

lbWlloo<IIPo||ooell(^^)+ll--*. (49) 

Let now K G N, K < k, i = (n, . . . , ^at) G and x G be such that YljLi H = K, 
\di-, . ..di^p{x,t)\ = \\di^ . .■di^p{t)\\^, and we define 

ll^'^'^llp = supIII^o-pII^p , \a\ < K}, l<p<oo,K <K, 

where a = (cti, . . . , ctat) G N^. W.l.o.g. we suppose that di^ . . . dij^p{x, t) > (to change the 
sign of this term, one just needs to replace the element ei of the basis of by — ei), and 
then, 

dtdi^ . . . di^p{x, t)=di^... di^Vx ■ {p{VW * p)){x, t) 

K 



E • (5ac(VW^*p)(x,t)) 



fc=0 \a\=k 
K 



+ E E(^-'^(^'^))(^-^(^^*'°)(^'*))• 



fc=o \(j\=k 

Here, a < i denotes aj < ij ioi j = 1, N and a'^ = a — i. Using that the term k = K in 
the first sum is zero one obtains 

dtdi^ . . . dij^p{x, t) =Vp{x, t) ■ (di^VW) * {di,, . . . di^p{x, t)) 
+ p{x, t){AW * {di, . . . di^p{x, t)){x, t) 

+ {d,, . . .di^p{x,t)){Aw * p){x,t) 

K-l 

+ E E (5aVp(x,t))(VW*a,cp)(x,t) 

fc=l \a\=k 
a<i 

K-l 

+ E E {dMx,m^W * Vd,cp){x,t), (50) 

fc=l \a\=k 
cr<i 
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and then, we get the estimate: 

dtdi,...di^p{x,t) < \\yp\\^\\di,VW\\j^i l|9*2 •••^iiv/'lloo 
+ \\di,...di^p\\^ IIAVFll^i IIpII^ 

K-l 

+ E E \\d^^p\\oo\\^W\\L^\\d^cp\\ 



k=l \a-\=k 
(7<i 

k=l \a\=k 



ivw^ll^i 



|V9o-cp|| 



< IIpIL II AVFlLi ll^i, . . . di,p\\^ + lia., VT^II^i ||9,, . . . d,,p\\^ \\V p\l 
+ Ck\WW\\^, 



W'^-^pf Wd'^pW 

I 1 1 oo M M oo 



An induction scheme on K initiahzed by ( 49 ) then provides the following exponential control 
for K G {0, 

Wd'^pW^ < Ci,Ke''^-'<' (51) 

where Ci^k,C2,k only depend on K, ||Z)^Ty||^i, ||(AVF)+||^oo, and ||/Oo|lw«'°°- Coming back 
to (50) the following estimate on the time derivative follows: 

d 



dt 



D^p 



(52) 



Finally, taking the derivative in t on ([T]) we get 

dfp ={VW * dtp) • {VW *p) + V^p • {VW * dtp) 
+ PiAW * dtp) + dtpiAW * p), 
from which it is easy to derive 



dt^' 



(53) 



Step 2: Construction of a solution trough an approximation problem. Let be a smooth 



are uni- 



approximation of W, that is G W^'°°(M^) such that ||D^W^||^i, ||(AM/^)+| 
formly bounded, and: 

VW'-^VW in L\R^). 

Thanks to 03111], there exists a classical solution p^ £ w}^^ {M+,W^'°°{M^)) with initial 
data po for each regular interaction potential . 

The estimate ( [sT] ) provides a uniform bound on ||/0^||oo! II^P^Iloo ll-^^'^^lloo' Since 
K > 2 then (52) implies that ^p'^ and (^Vp^) are uniformly bounded for e > 0. Applying 
the Ascoli-Arzela theorem, due to ( [sT] ), ( [52] ) and ([53]) there exist limits for p"^, dtp'^ and V^p^ 
(where we have written e instead of e^) on C([0,T] x B) for any compact subset B C 
and moreover the limits denoted by p, dtp and V^p belong to C([0, T]; W^'°°(M^)). For the 
velocity field we have that 

\v^x,t)-v{x,t)\< [ \{VW' -VW){x-y)\p'{y)dy+ [ \VW{x - y)\\p'{y) - p{y)\ dy 

= {I)e + {II)e- 
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For (/)£ one observes that 

|(/)^| < ||VW/^-VW^||^i ^0 as £^0, 

and for {II)e one has that that 

\VW{x - y)\ \p%y) - p{y)\ < C \VW{x - y)\ E L\R^). 

Then, by the dommated convergence theorem, {II)e — )• as e — )• and, as a consequence, 
v^{x, t) converges pointwise to v{x, t) in x [0, T] for ah T > 0. The same reasoning is used 
to prove that V^^ • — )■ V^; • w as e — 0. Thus, the regularized equation 

= V/>" • {VW * p') + (AVF= * p')p'') 

passes to the hmit and p E CH[0,T] x M^) n W^'~([0, T], W^'°°(M^)). 



The propagation of the regularity follows from estimates (^51| and (52). The proof of 



uniqueness follows from |23| [TT] . □ 

Remark 8. Under the assumptions on W in the previous theorem we have that p is Lipschitz 
continuous both in space and time, and then the characteristics are well defined: 

^^Xt = -{VW*p){Xt,t), 

and the solution p is given by 

pix,t) = po{Xf')detiDX^'). 

Remark 9.IfW and po are radially symmetric, then one can easily check that the problem 
is invariant through rotations around the origin. The uniqueness result then shows that the 
solution p is radially symmetric at all times. We also point out that if the solution is compactly 
supported then it remains of compact support for all times. 

6. The example of power law repulsive-attractive potentials 

The aim of this section is to show an example of how to apply the general instability and 
stability theory in the case of power law repulsive-attractive potentials: 

I la 11^ 

Wix) = ^_^ 2-N <b<a. (54) 
a b 

The condition b < a ensures that the potential is repulsive in the short range and attractive 
in the long range. One can easily check that for these type of potentials AW E L^^(M^). 
The condition 2 — N < b ensures that the potential is in Wlf^{R^) for some I < q < oo. 
Using algebraic computations, involving the Beta function, we give the conditions that the 
powers a and b should satisfy in order to apply the stability and instability theory, and we 
construct the bifurcation diagram for these powers. The main results of this section are the 
following: 

Theorem 8 (Global existence of solutions for repulsive-attractive potentials). Given W by 
(54). Assume po E yV'^'°°{M.^) is compactly supported and radially symmetric. Then there 
exists a global in time classical solution for Furthermore, the solution is compactly 

supported and confined in a large ball for all times. 
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Theorem 9 (Sharp radial stabihty-instabiUty for spherical shells). Assume that W is a power 



law potential as in (54). Then, there exists a unique Rab > given by 

1 //3(^^±f^ 




"2 v/^c^,^), 

such that 5r^^ is stationary solution to ([T]). Moreover, the following properties hold: 

(i) If2-N<b<3-N then oj £ C{Rl) n C\RI \ V) and for all {R, R) £ V we have 

lim (r, 77) = +00. 

(r,J7)^(H,i?) dr 

(ii) Ifb£ (s-N, 3"-^°-^°+^3^-^' ) then Lo is Ci(R^) and 

diUj{Rab,Rab) > 0. 
(Hi) If be |^ 3a-Ara-10+7jV-Af2 ^^^ ^f^^^ ^ C\RI) and 

diUj{Rab, Rab) < and {diOJ + d2Uj){Rab, Rab) < 0. 



As a consequence, if b £ {^2 — N, ^ ^ j then 5^^^ is unstable in the sense of 

Theorem jsj and if b £ ^ 3a-Afa-40+7JV-jV then Sr^^^ is stable in the sense o/ Theorem joj 

Remark 10. Note that indeed for 3 — N < a we have 

3a- Na- 10 + 7N - iV^ 

3 - iV < — < a. 

a + N -3 

Remark 11. In |39j the authors study the dynamic of a curve evolving in R? according to 
the aggregation equation. They perform a linear stability analysis of the spherical shell steady 
state. They consider not only radially symmetric perturbations but also perturbation which 
break the symmetry of the spherical shell. The mode m = oo corresponds to a perturbation 
which preserve the symmetry of the spherical shell. Using a computation involving the Gamma 
function, they show that the mode m = oo is stable if and only if (a — l){b — 1) > 1. In order 
to prove (ii) we will perform similar type of computations involving special functions. Note 
that, when N = 2 

3a- Na-lO + lN - N"^ _ a 
a + N -3 ~ a- 1 

so is equivalent to (6 — 1) (a — 1) > 1 and we recover the condition derived in |39j . 

As a summary of all the stability and instability results for a Sr^^^ stationary states for 
power law potentials we show the bifurcation diagram in Figure [l} For powers inside the 
region between b = 2 — N and the curve one has instability of the 5Il^^^ . In the region above 
the curve one has stability of the 5r^^^ . 
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6.1. Proof of Theorem [8j In order to prove the global existence theorem we need to in- 
troduce some notations. For potentials defined by (54), the kernel uj{r,ri) defined in ^ 
becomes 

a;(r, rj) = r^~^iljiy{rj/r) — r°'~^ipa{r]/r) (55) 

Ms) = i / Mv) (56) 

f^N JdB{o,s) \ei - sy\ 

The properties of the function V'a(-s) that we need are summarized in the following lemma 
and can be found in [30l. 



Lemma 7 (Properties of the function ipais))- The function ipa is continuous with tl'aiO) = 1 

+a- 
N 



and lim.^oo s^-^M^) " ^^^^ 



The main difficulty we have to cope with is the growth at infinity of the attractive part 
which restricts the range of direct application of Theorem [7] 



Proof of Theorem^ Due to translational invariance we can assume without loss of generality 
that the center of mass is located at zero. We write W as W{x) = Wfi{x) + Wa{x) where 

Wa{x) = is the attractive part and Wii{x) = —-^ is the repulsive part. In addition, since 
W is radially symmetric, that is W^(a;) = A:(|x|) then we define k{r) := kji{r) + fcyi(r). Finally, 
we write u;a('^, = r"'~^ipa{il/''^) ^-nd a;jj(r, r/) = r^~^il)i){r)/r) withuj{r,ri) = uiji{r,ri)—ujA{r,r]). 

Step 1: A priori estimates on the support of p. Suppose that p is a smooth radially sym- 
metric solution for the equation ([T]) with compactly supported initial data po ^ W^'°°(M'^). 
Since the solutions belong to W^^^J^ (M-(_ , >V^'°°(M^)) and are compactly supported, the veloc- 
ity field V is Lipschitz continuous in time and space and the characteristics are well defined. 
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Thus V generates a flow map r(t,ro), t £ [0, T], tq G M+: 

r(0,ro) = ro. 

Let us define by r2{t) the characteristic curve starting at point r2(0) = max{supp(/io)}- Then, 
—r2{t)=v{t,r2{t)) = J^ uj{r2{t),'n)dfit{v) = ^{'r2{t),ri)dfit{r}) 

r2{tr-^Mr]/r2{t))dfit{v)+ r2{t)''-^Mv/r2{t))dflt{ri), (57) 



where we have used the expression of uj given in ( |55[). Here fit denotes the measure with 
density pt- Using the properties of il^a in Lemma [7] in (57) we obtain the following inequality: 

"^Mt) < KbT2{tf-^ - Kar2{tr-\ (58) 



where 



dt 

Ka = 1 and Kb = V'b(l), if 2 < 6 < a, 

Ka = l and Kb = l, if 2 - iV < 6 < 2 < a, 

Ka = 'ilJa{l) and Kb = l, if 2 - iV < 6 < a < 2. 

Defining Rab '■= (^"^^ " and rewriting (58) as 

j^r2{t)<r2{tr-\KbT2{tf-''-Ka) 

one realizes that r2{t) < R := max(r2(0), Rab) which proves that the supp(/i() is bounded and 
contained in .6(0, R) for all times. 

Step 2: Global existence. Given < e < 1, consider Xeii^) a C°^(0,oo) decreasing function 
with < e < 1 such that Xe('^) = 1 if < r < 1/e and Xe{f) = if r > 1 + 1/e. Define 
/^(r) := xe(r) • k'^{r) £ L^0,oo), k%{r) := Jq f%s)ds and k%r) = kR{r) + k%{r). Now, the 
potential W^{x) = k'^{\x\) satisfies the hypotheses of Theorem [7] so we have existence and 
uniqueness of classical solution to ([5])-(|6]) in [0,T] with initial data po. We denote by pf 
the measure with density pf. 

Consider r2 = max{supp(/io)} and the characteristic curve r^f) starting at point r2 = 
r2(0). Computing the derivative with respect to time, one has 

-rl{t) = uj'{rl{t),r,)dp'M = u;%r'^{t),r])dpt{rj) 

u%{rm,v)dpt{v)+ / rm'-'Mv/rm)dpKv) 
Jo 

< Cbrl{tf~\ 

where we have split the kernel uj'^ into its attractive and repulsive parts, and we have used 
that (j^ > 0. The constant Cb depends on h. The last inequality leads us to 

rm<'y{t) ■.= {Cb{2-b)t + rl-')^\ (59) 
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which says that the solution exists at least up to time T* := ^ min{r, T^}. where 

I +00 if 6 < 2. 



In addition, ( |59[ ) gives us a uniform estimate for the support of up to time T*. Notice that 
for aU f < T* and e > such that 2a{T*)£ < 1, then VW^ * Pt = VVF*pf for aU x £ supp(/of) 
and all t G [0,T*]. As a consequence u given by (55) and u'^ associated to by ^ are 
equal in the set {(r, 77) | {r,7]) £ supp(/9f)^} for all t < T*. Therefore, we can write: 



dt 



a;^(r|(t),r?)d/if(r?) + 



rm 



rm 



WA(ri(t),r?)fi/if(r?) + 



rm 



rm 



rm^-'Mv/rmWliv) + 



rrm 



and we can use the a priori estimates developed in Step 1 . Then, for all t <T* we can conclude 
that r|(t) < R. Now, let us take £ such that 2£R < 1. Therefore VW^ * pf = VVF * pf . For 
all t < T* in the support of pf. By uniqueness pf =: pt for all 2eR < 1 and it is a classical 
solution to ([1]) with potential W. Summarizing, we have shown the existence of solution in 
the time interval [0,T*] with r2(T*) < R. Now, we can extend and repeat this argument for 



a time step At := ^ min (l,Tfe), where Ti, 



if 6 > 2 or Tf, = +00 if 6 < 2, obtaining 



a solution up to time T* + At such that r2{t) < R for all t G [0, T* + At]. Since At is 
independent of the initial data and e > 0, then we can extend the solution for all times. 

Finally, the a priori estimates on the support of p show that the support of the solution 
remains compact for all times. □ 



6.2. Proof of the Theorem [oj It is first convenient to rewrite (56) as: 

aN^i r cos 0) (sin 0)^-2 



UN 



2s COS / 



2-a 



de. 



(60) 



/O (1 + s2 

We recall that i^(r, r/) is the velocity at r generated by dB{0,r]). So a 6r with > is a 
steady state if and only if uj{R, R) = 0, i.e. 



R — R. 



ab 



V'a(l) 



1 

a-b 



(61) 



where we have used (55). 



Proof of the Theorem [9[ The point (ji]) is a direct consequence of Lemma [Sj Let us prove ([n]) . 
From Lemma [sj see also |30| it is clear that uj £ C^(M^) and we have 



dr 



{Rab, Rab) = Kf (b - 1)V6(1) - 4^'bil] 



R 



.a-2 
ab 



i)V^,(i)-<(i) 



(62) 



After some algebra, one easily get from (61) and ([62j) that ^ {Rab, Rab) > is equivalent to 



V'a(l) 



< 6 



mi 



(63) 
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Both 'i/'a(l) and i^'ai^) can be expressed in terms of the Beta function. Recall that one of the 
expressions of the Beta function is: 

»7r/2 



Jo 



We first compute '0a(l)- Using (60): 
(Tat , /I \ f^i^~ cos 



Ml)= r^^v^^(sin<-2d0 = 2^ r{l-cos0r/Hsin9)^-'de 
Jo A[L,U) Jq 

2^ r(2sin2-r/2(2cos-sin-)^-2<i^ 
Jo 2 2 2 



/■TT / /i\ a+Af-2 / a\ N-2 



where we have used the fact that A{1,9) = y^2(l — cos 0) and the identities 1 — cos 6 = 2 sin 
and sin = 2 cos ^ sin ^ • Similarly we compute 



2 e 

2 



a^_i(a-2)(a + iV-2)"«^ ^ 7o ^(1,^)"-'^ 7o (2(1- cos 0)) 



2 



- (2 COS I sin f)^ 



,2 



(2(2sin^|))- 

N / Q\N+a-4: 

COS - 1 I sin - I d9 



•,N+a-A 



/ cos - sm - 
Jo \ 2; V 2; 



Note that since a + — 3 > the Beta function is well defined. If we compute the quotient 
we obtain: 

V^;(l) _ 1 (a - 2)(a + AT - 2) /3 (^, ^) 



V'a(l) 2 iV-1 /3(«+A^,A^)- 

r(^)r(t) 
r(2+t) 



At this point, we remind that /3{z,t) = ^^fz+^t) ■ With this expression, the quotient can be 
simplified as 

V.^(l) _ 1 (g - 2)(a + AT - 2) r (^) r (^) 

Va(i) 2 AT-i r(^^±f^)r(^)' 

and if we use that T{z+1) = zT{z) then, the gamma quotient can be reduced to (A'' — l)/(a + A^ 
and then we obtain 

V.;(l) ^ l (a-2)(a + A^-2) 
V'a(l) 2 a + A^-3 



Plugging the above expression into (63) and doing some algebra we deduce 

{a + N- 3)52 + [N^ -7N + 10- a'^)b - {N^ -7N + 10)a - (A^ - 3)0^ > 0. 
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The roots of the quadratic form are 

, , 3a- Na-10 + 7N - 
b = a and b = 



a + N -3 

which gives To prove (iii) one just need to replace the < sign by a > in (63) to get the 

dm ~~ 

dr] 



first condition. For the second condition of stabihty, we can easily compute ^ from (55) and 



(|62|) to get 

+ ^) (^'^'^ ^«^) = ^afib - ma) - Kf{a - l)Ml) 



Now, using the definition of Rab in ( |61| ), we finally obtain 

' duj duj\ , 

- I {Rab, Rab) = b-a<0. 



dr drj J 

The stated instability and stability are direct applications of Theorem [3] and Theorem [6] 
respectively. □ 

Remark 12. We want to point out that the general theory developed in the previous sections 
is still working in dimension N = 1 for even solutions which correspond to the radially 
symmetric solutions in higher dimensions. In the case N = 1, and for even solutions, the 
function V'a; corresponding to W{x) = reads 

Ms) = ^ [(1 - s) |1 - sl"-^ + (1 + s) |1 + s| 

One can easily check that the properties of the function il^a and uj for W(x) = '—^ '—^ in 

N = 1 are the same as in Lemma [t] and in Theorem [oj The radius is Rab = | whatever the 
powers are, see (61). Theorem [9] app/ies.' if b ^ (1)2) then we are in the instability case " 



and if b G [2, a) we are in the stability case. The curve which separates the instability and 
stability regions in Figure^ degenerates and becomes the line 6 = 2. 

Moreover, in [331 [32] the authors proved the existence of weak solutions and convergence, 

up to extractions of subsequences, of p{-,t) for potentials like W{x) = b G (0,1]. 

They also showed numerical simulations supporting the conjecture that the stationary state 
I ((^2;=_i/2 + 3^=1/2) 'is unstable. Note that our Theorem [oj only applies in N = 1 for 1 < 
b < a. Since global existence was proven in [33\ I32j . then our instability result also applies 
in those cases. Summarizing, we can include a = 2, 6 G (0, 1] in the instability regions using 
Theorem [H 

7. Numerical results 

In this section, we illustrate the previous results and get some further conjectures for the 
instability cases. Our numerical code is based on the inverse distribution function in radial 



coordinates. As it was reminded in (31), the equation for the inverse distribution function 
reads 

^{t,0= [ oo{v{t,0,v{t,i))dl (64) 







A solution of ([s]) converges to a Dirac mass if and only if its pseudo inverse distribution 
becomes flat. 
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Numerical codes based on (64) are then more stable when dealing with mass concentration. 



We will then use a backward Euler scheme in time coupled to a composite Simpson rule to 
approximate the integral term, and solve the resulting nonlinear system by the Newton- 
Raphson algorithm. Let us remark that the convergence of the semi-discrete backward Euler 
scheme is equivalent to the convergence of the JKO variational scheme for ([T]) (see |13| I22j). 
The convergence of the semi-discrete backward Euler scheme is therefore known under suitable 
conditions on the interaction potential, see [18j for details. All simulations are done for N = 2. 



In this case, (64) reduces 



Test case, total concentration at the origin: W{x) 
to + = 0. To test our scheme, we use this attractive potential for which the 

solution converges exponentially fast to a total concentration at zero, that is to = 0. See 
Figure [2] 




0.2 



Figure 2. Evolution of C ^ ^it,0 for W{x) = ^ 
tion at 0. 



towards total concentra- 



Stability Case for the Spherical Shell: W{x) = In this case, we have an 

repulsive-attractive power law potential with powers in the stability region of Figure [T] We 
thus expect that the mass will concentrate towards a spherical shell, thanks to the results of 



Theorem [9j The radius of the spherical shell can be computed using (61): 

„ _ fM^)V^ _ V3 

"""'-[MI)) " 3 • 
For b >2 and both a and b integers, one can compute explicitly the expression for the velocity 
field uj{r,r]), which is a polynomial function, in our case uj{r,rj) = —r^ — 2rrf + r. The 
evolution of ip is shown in Figure [sj In Figure [s] we also plot the velocity field r i— )• a;(r, Rab)- 
Notice that r i— )• u!{r,Rab) satisfies the conditions of Theorems [s] and [6j oj{Rab,Rab) = 0, 
diUj{Rab, Rab) < 0, sign{u; {r, Rab)) = sign{Rab - r), diUj{0,Rab) > 0. 

Instability Case for the Spherical Shell: W{x) = In this case, the powers 

are in the instability region of Figure [T| below the curve b = Then, due to the results 
in Theorem [91 a spherical shell is unstable. One can notice on Figure 111 that the function 

— I |2 — 

r I— >• uj{r,Rab) associated to the potential W{x) = |x| satisfies diuj{Rab,Rab) > 0, so 

that the instability condition of Theorem [3] is indeed satisfied. 
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n.2 0.4 O.fi 0.8 1 0.1 0.2 0.3 0.4 0.5 0.8 0.7 0.8 

Figure 3. Case W{x) = - i^. Left: Evolution of ^ towards 

the uniform distribution on the sphere of radius Rab = Right: Velocity 
field r i— )• uj(r, ^) with the vertical line pointing out Rab- 




Figure 4. Case W{x) = '-^ Left: Evolution of ^ i-> ip{t,(,) towards 

a stationary profile, possibly an integrable function. Right: Velocity field 
r I—)- uj{r,Rab) with the vertical line pointing out Rab ~ 0.6366. 



Figure |4] shows that the solution seems to converge to some stationary state which does 
not have any singular part, i.e., possibly an integrable function. Numerically, this behavior 
appears for any powers a, b in the instability region of Figure [T] We conjecture that in this 
region there exists integrable radial stationary states which are locally stable under radial 
perturbations. This has already been proved in the particular case oi b = 2 — N and a > 2 
in |34|- Some numerical simulations using particle systems done in [39] however suggest that 
these stationary states might be unstable for non radial perturbations. 

Energy dissipation. We remind that the energy functional is given by 
E[p] it)= [[ Wix - y)pit, x)pit, y) dy dx . 
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Using the polar change of coordinates x = ra and y = sa and using the radial symmetry of 
p{t, •), this energy writes: 

E[p]{t) = -^ll I W{ra - sei)p{t,r)p{t,s)dadsdr . 



2c^Af J Jr2^ J 83(0,1) 




10 12 14 



Figure 5. Energy decay in logarithmic scale for the regular repulsive- 
attractive potential, case a = 4 and 6 = 2 (solid line) and for the singular 
repulsive-attractive potential, case a = 2 and 6=1 (dashed line). Note that 
Emin is the numerical limit of the energy as t — )• oo. 

A formal calculation implies that the derivative w.r.t. time of the energy is negative and 
given by 



d^ 
dt 



E[p]{t) 



p{t, r)v{t, r)^ dr . 



the energy should then decrease in time. Using radially symmetric coordinates, the energy 
functional for the inverse distribution function is given by 



Em) 



1 



2a N 



1 fi 



Jo JdB{0,l) 



W{Lp{t, - if{t, ^)ei)dadidi . 



(65) 



We have computed the energy using the formula (65) to check numerically, in each case, that 



the energy decreases. In Figure [5] we observe the exponential decay of the energy for the two 
numerical examples presented above for repulsive-attractive potentials. 



8. Appendix 

Let us start by some differential geometry facts. For the sake of clarity, we first define the 
type of hypersurfaces we will work with. 

Definition 6. A4 C is a hypersurface (manifold of dimension N — 1) if for any x ^ A4 
there exists a chart {U,ip), i.e., a pair of an open connected set and a C"^ diffeomorphism 
if : U — > M^, with x e U C R'^ such that ip{x) = and y e M D U if and only if 
^{y) G {0} X 



We will need some technical result from differential geometry in order to deal with the 
regularity of the function w in ([2| and its generalizations to any compact hypersurface. Note 
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first that ii A4 C is a hyperplane then M n dB{x,r) is a — 2 dimensional sphere of 
radius (r^ — dist(j;, A^)^)^^ and therefore its surface area is 

JV-2 

\M n dB{x,r)\-^N-2 = cTN-iir"^ - dist(x, A^)^)^^ 

where H.'^ is the d- dimensional Hausdorff measure, and we remind that ittv-i is the surface 
area of the unit sphere in M^~^. 

The following result is a classical consequence of uniform graphs lemmas in differential 
geometry. They state that a compact regular hypersurface can be covered by graphs with 
bounds on their derivatives depending only on the uniform bound of the second fundamental 
form. We refer to [SUI Lemma 4.1.1]. This allows to show that the volume elements locally 
converge to those of a hyperplane in a uniform manner. 

Lemma 8. Let A4 C be a compact hypersurface of dimension N — 1 immersed in 
M^. Then there exist small enough ro > and constants C,C > depending on the global 
bound of the second fundamental form of Ai such that for all < r < ro, and all x G 
with dist{x,M.) < ro 

N-2 N-2 

(7(r2 -dist(x,Ai)2)+2 <\MndB{x,r)\-^N-2 <C{r^ -dist{x,Mf)+^ . (66) 

Remark 13. Let us note that the previous Lemma is trivial in the case of M = dB{0,r]) 
for any rj > since the intersection of two {N — 1)- dimensional spheres of different radius is 
always a {N — 2) -dimensional sphere lying on a hyperplane. In fact, we can easily compute 
that if two spheres dB{0,ri) and dB{x,r) intersect, that is \\x\ — ?]{ < r, then 

\dB{0, T]) n dB{x, r)\^N-2 = fj^v-irf 

where ri = ri(ry, r, dist(x, 5i?(0, r/))) is the radius of the intersection, which is computable: 



ri 




2|x|r/ 




13 



can then be taken 



as r — dist(x, 9-6(0, r/)) — t- 0. The constants r^, C, and C of Lemma 
uniform for variations of the radius in bounded intervals, i.e., for < r/i < rj < r]2. 

We now can deal with the continuity of the velocity fields generated by probability densities 
concentrated on manifolds. Recall that M+ = (0,+oo). 

Lemma 9. Let Ai C be a compact hypersurface, fi a probability distribution such that 
j2 = 4>6m; where (j) G L°°{A4), and g E C(]R^/{0}) a radially symmetric function which is 
locally integrable on hypersurfaces. Then, the function 



^(^) = / 9{x-y)dfi{y) 

is continuous in x ^ R-^. Moreover, the same results hold while replacing g{x) by a non- 
radially symmetric function G S C(M'^/{0}) such that |G(x)| < |5(x)|, where g satisfies the 
properties above. 

Proof. It is straightforward to check that v[x) is continuous for all x ^ M. Let x ^ M and 
let ro be given by Lemma Isl For < e < ro, let Xe £ C°°(M+) be a cut-off function, such 
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that Xe = 1 on [0,e/2], and Xe = on [e,oo). The function v can then be written as 

v{x)=\ g{x - y)xe{x - y) dfl{y) + g{x - y)[l - Xsix - y)] dp,{y) 
JM JM 
■.= vl{x) + vl{x). (67) 

It is clear that continuous on x G M'^, since g is continuous away from the origin and 
supp [jl) = is compact. Moreover, given the set C/ = {x G : dist(x, A^) < ro}, we can 
estimate for all x G f7 

\vl{x)\< I \g{x - y)\xe{x - y)dfl{y) < \g{x - y)\ dp.{y) 

Jr^ JB{x,e) 

<U\\l^{M) [ \9{r)\\{y(^M;\y-x\=r}\^N-2dr<C\\(P\\L^(^M) [ \9ir)\ r^''^ dr 
Jo Jo 



where (66) is used. Moreover, by construction wf(x) = for all x ^ C/ for e < tq. Therefore, 



due to the integrability over hypersurfaces of g, then 

lim \\vf\\Lc^(^f.N) =0. 

This is enough to show the continuity of u on A^: for any 6 > 0, there exists e > such that 
11''^! IIl°°(r'V) ^ |- Since is continuous, there exists k > such that I^^IC^) ~ ^1(^)1 — I 



|x — x| < K. Then, \v{x) — v{x)\ < 5 if |x — x| < ac. The last part of the proof is an adaptation 
of the previous arguments since the integral inside the norm is less or equal than vf'' . □ 

Now, we want to obtain the continuity with respect to the hypersurface for the velocity 
fields associated to measures concentrated on them. We restrict to the case of spheres since 
we only need this particular case. The proof uses the transport distance doo- We remind the 
reader that it is introduced in Section 3. 

Lemma 10. Let := dB{0,rj) and Jlrj = ip'n^Mr, probability measures such that c/)^ G 
L°^{M.^) with < T]. Let g G C"^(M'^\{0}) be a radially symmetric function which is locally 
integrable on hypersurfaces. If the functions are uniformly bounded in tj and dooifir], fifi) — ^ 
as rj — fj ^ 0, then 

v{x,r]) = / g{x - y) dfiniy) 

JdB{0,rj) 

is continuous in x M4-. Moreover, the same result holds while replacing g{x) by a non- 
radially symmetric function G G C^(M^\{0}) such that |G(x)| < \g{x)\ with the properties 
above. 

Proof. Lemma [9] implies directly the continuity with respect to x for all fixed 77. Using the 
Remark 13 and the proof of Lemma [9j it can be easily checked that this continuity in x 



is uniform in r/. Indeed \vW can be made small uniformly in rj and, due to the estimate 
|Vu|(x)| < supg5(^ l^bll ~ ^2 is coutinuous uniformly in rj. Therefore, we only need 

to show the continuity in 77 of for a fixed x G M^. 



As in the proof of Lemma^ let rg be as obtained in Remark 13 uniform in < 771 < rj < r]2. 



We choose again < e < tq and Xe £ C°°(M+) a cut-off function, such that Xe = 1 on [0, e/2]. 



and Xe = on [e,oo). We can write v{x,r]) = vf{x,r]) + V2{x,rj) analogously to (67). As in 
Lemma [9] using the properties of g and the uniformity in Remark (13), we can easily show 
that 

lim ||uf(-,r/)||^oo(KiV) =0. 
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uniformly in < r/i < r] < r]2- Therefore, for any 5 > 0, there exists tq > e > such that 
!!''-'!(■) ^) IIl°°(m^) ^ f uniformly in < ?7i < < 772- Now, we estimate 

\v{x,'q)-v{x,i))\ < \\vl{-,r])\\L^^f,N) + ||^^f (•, ??)||L-.(iRiV) 



+ 



gix -y)[l- x%\x - y\)] difir, - i^f,){y) 



< 2 + - X^)]||l°°(A) doo{flr^,flfj), 

where A is the convex hull of the set {x} — (supp firj) U (supp ftfj). Notice that the set A is 
uniformly bounded in r] and fj. 

This estimate shows the continuity in r] since dooip-r], fifj) — ?• as r/ — )• 77, and thus, the last 
term is bounded by (5/2 provided that ry is close enough to fj. Again, the final part of this 
Lemma is a small variation of the previous arguments. □ 

Finally, we complete the results by showing that if the function is not locally integrable on 
hypersurfaces then the velocity field is not bounded. 

Lemma 11. Let Ai^i '■= dB(0,rj) and firj = 4>r]SM,^ be probability measures such that (prjix) > 
> for all rji < rj < ri2. Let g E C(M^\{0}) be a nonnegative radially symmetric function 
which is not locally integrable on hypersurfaces. Then For all M > there exists 6 > such 
that 

dist(x, Mr]) < S =^ / gix — y) djir^iu) > M for all x G and for all 7?i < ?? < f?2- 
Proof. Using Lemma [s] and Remark 13, for x G with dist(j;, A^r;) < '''o, we get 



/ 



g{x - y) djlriiy) > j g{x - y) dflr,iy) > 4>o gir)\{y € Mr,; \y - x\ = r}\y^N-2 dr 

J \x—y\<ro Jo 

f-ro jv-2 

>4oC g{r){r^-dist{x,Mr,f)+' dr. 
Jo 

Since Jq g{r)r^~^ dr = +00 and g is continuous and nonnegative on (0, 1], we deduce that 

rro jv-2 
lim / g(r) (r'^ — dist(x, Mri)'^):^ dr = +00 , 

dist(a;,M^)^0 Jq 

by the monotone convergence theorem, which conclude the proof. □ 
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